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Метод выделения областей на изображениях 
с использованием фрактальных и текстурных 
характеристик снимков высокого разрешения 


В статье предлагается метод обработки изображений, основанный на текстурных и фрактальных 
характеристиках изображений. Приведены результаты совместной сегментации с использованием 
разработанного метода при обработке аэрофотоснимков сельскохозяйственных полей и схема программно- 
аппаратной системы принятия решения, использующая разработанный алгоритм обработки. 


Введение 


Дистанционные методы мониторинга сельскохозяйственных полей дают воз- 
можность оперативно выявить участки полей, пораженные болезнью. Выявление за- 
болевания на ранних стадиях развития значительно сокращает затраты труда и 
средств. Находят развитие два основных подхода к решению задачи выявления 
пораженных участков — спектрометрический и оптический [1], [2]. Спектральный 
подход позволяет определять многие заболевания на ранних стадиях развития. Однако, 
несмотря на это, оптические методы также существенно развиваются, так как их 
использование позволяет увеличить качество идентификации. 

Объектом исследования являются цветные изображения сельскохозяйственных 
полей (рис. 1). 

Цель работы состоит в разработке эффективного метода обработки полутоновых 
и цветных изображений растительных покровов, полученных с помощью цифровой 
съемки высокого разрешения, а также их реализации в виде программных средств 
для систем компьютерного зрения. Суть метода состоит в разработке и использовании 
фрактальных и текстурных характеристик с требуемыми свойствами для построения 
пространства признаков и выделения объектов на цветных изображениях сельскохо- 
зяйственных полей. 


Текстуры 


Неоднородность либо повторяемость мелких фрагментов называется текстурой 
цифрового изображения. Характерным признаком текстуры является однородность 
на уровне окрестностей или локальном уровне, т.е. на уровне групп смежных пикселей 
с различными значениями яркости [3]. 

Текстурные изображения принято делить на два класса: стохастические (или 
случайные) и периодические (или структурные) текстуры. На основании признаков 
базовых элементов они подразделяются на мелкозернистые, крупнозернистые, гладкие, 
гранулированные и холмистые. 
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На основании степени взаимодействия базовых элементов подразделяются на 
сильные (взаимодействие подчиняется некоторому правилу) и слабые (взаимодействие 
имеет случайный характер). 


— 
— 


‚8 


у я Ш мч = < 


100 метров —_ 50 метров 


15 метров 5 метров 


Рисунок 1 — Примеры исходных аэрофотоснимков 


В результате анализа основных текстурных характеристик для использования 
выделены следующие: АЗМ (Апецаг Зесопд-Мотеп! еааге) — мера однородности 
изображения; Сопёгаз{ — мера количества локальных вариаций на изображении и 
Епгору — мера беспорядка в яркостях пикселей изображения: 


двм- 55268, (1) 


1=1 7=1 


где № — количество отдельных серых уровней на изображении; Р - матрица 
относительных частот; К — количество пар пикселей, 
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Выбор указанных характеристик основан на результатах проведенного визуаль- 
ного экспертного анализа и опирается на наличие видимых различий между больны- 
ми и здоровыми участками полей. Пример результата вычисления этих текстурных 
характеристик приведен на рис. 2. 


Рисунок 2 — Результат вычисления текстурных характеристик 
(в данном случае — АЗМ) для аэрофотоснимка поля, 
выполненного с высоты 50 метров 


Фракталы 


Суть вычисления фрактальной сигнатуры состоит в том, что квантованные зна- 
чения интенсивности двумерного сигнала должны располагаться между двумя функ- 
циями, называемыми верхней и нижней поверхностями. Верхняя поверхность ИЦ 
содержит множество точек, значения которых всегда, по крайней мере, на один 
квант интенсивности превышают интенсивность входного сигнала. Нижняя поверх- 
ность / имеет значения точек, которые всегда ниже, по крайней мере, на один квант 
интенсивности входного изображения. 

Верхняя и нижняя поверхности при нулевой шкале масштабов определяются как: 


О (6, 1,0) = Ь 1,0) = #07}; (4) 
где 2(71, /) — входное изображение. В общем случае: 


ОС, = +1) = шах 1-5) + пах (вт, 2) (5) 
Вьет 0,6) рат, ©) 
п = { (т) [Ч [(®,т), (6 151; (7) 


где 4- функция расстояния. 

Сконструированное покрытие, образованное двумя указанными функциями, 
имеет толщину 2 = . Для двумерного сигнала площадь поверхности есть объем, занятый 
покрытием и деленный на величину 2=. Площадь «поверхности» интенсивности 
А(=) в приделах окна наблюдения А рассчитывают вычитанием точки за точкой 
нижней «поверхности» из верхней с дальнейшим суммированием по всему окну: 


У Ча, =)-МЬЛ=) 


Ав) = в 


2= 2= (8) 
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Фрактальную размерность определяют по наклону [02 А(= ) как функцию /09=. 
Пример фрактальной сигнатуры изображен на рис. 3. 


Рисунок 3 — Фрактальная сигнатура участка изображения 


Фрактальная размерность 0(1,/) при нахождении пикселя (1,]) на всех 
масштабах оценивается как взвешенная сумма локальных фрактальных размер- 
ностей Р. (1, ]) в виде: 


Усе 
р) = Ус: (9) 
при | 
@= 105=-05(=-—1. (10) 
1052 
к ООО ЮВ ВО: (1) 


105 &-105(=-1) 


Вследствие свойства, изложенного в [4], величина Р. (1, /) является делением 
А(Ь Е ) на А(1, , 5 - р: 


А. ),=) _ Кё” -[ в” (12) 
АС ,Е- К@Е-б” \=-1 
Логарифмируя, получаем 


102 = —105(= —1) 


Подставив в выражение (13) значения слагаемых (11) и (12) и сократив 
слагаемые, получим: 


105 А(1, },=) — 05 (т, },1) 


О(р, 7) = 
7 105 = -—1051 


(14) 
С помощью выражения (14) можно строить оценки фрактальных размерностей. 
Результаты работы алгоритма вычисления фрактальных сигнатур представлены 


на рис. 4. 
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Рисунок 4 — Фрактальные сигнатуры различных областей изображения 


Метод обработки 


Суть метода обработки состоит в совместной сегментации [5], используя их в 
качестве дополнительных информационных каналов, дополняющих имеющийся 
исходный снимок. 

Вычисление фрактальных сигнатур и текстурных характеристик изображений 
осуществляется для отдельных каналов с последующим их объединением с использова- 
нием коэффициентов, значения которых зависят от типа и состояния растительности. 

В качестве пространства признаков, на основании которых принимается реше- 
ние, используются взвешенные матрицы цветовых характеристик исходного изобра- 
жения, а также текстурные и фрактальные характеристики, вычисленные для каж- 
дого цветового канала исходного снимка. 

В качестве цветовых признаков используются определенные экспертом цвето- 
вые диапазоны, соответствующие здоровым и больным участкам полей. В табл. 1 при- 
ведены значения диапазонов цвета больных участков растения. 


Таблица 1 — Значения диапазонов цвета 


Тип сегмента Диапазон оттенка Диапазон насыщенности 
Здоровый (зеленый) 1,4; 3,14] [40; 200] 

Больной (желтый) [0,9; 1,4] [80; 200] 

Больной (буро-зеленый) 1,5; 1,8] [14; 55] 


Значение оттенка вычисляется по формуле 


о 2 о .Ь 
Ние = агс ап 3 ; (15) 
2 2 
+ а п). 2+ а п)-Ь 


а насыщенности — по формуле 
за! = тах(г, е,Б) — тш(г, 2,6), (16) 


где г, ©, Б — значения красной, зеленой и голубой компоненты цвета. 
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Таким образом, можно предложить следующий алгоритм обработки исходных 
изображений для составления карт особых областей (рис. 5). 

1. Обработка исходных изображений для получения дополнительных информа- 
ционных каналов, представляющих собой матрицы текстурных и фрактальных характе- 
ристик каждого цветового канала исходного изображения в отдельности. 

2. Выполнение совместной сегментации полученных матриц текстурных и фрак- 
тальных характеристик и цветовых каналов исходных изображений. 

3. Построение карты особых областей на основании поставленной задачи и резуль- 
татов проведенной совместной сегментации. 


Текстурные 
характеристики 
Фрактальные Принятие 
характеристики 
РАВЕ 9 Совместная решения Карта особых 
Е № 


сегментация областей 


Исходные снимки Зы 


Рисунок 5 — Схема алгоритма обработки 


Данный алгоритм предназначен для выполнения сегментации исходных двумерных 
данных, представляющих собой матрицы различных характеристик исходного изо- 
бражения (в нашем случае в качестве этих признаков используются цветовые ка- 
налы, а также текстурные и фрактальные характеристики). Таким образом, работа алго- 
ритма осуществляется в №-мерном пространстве признаков (где №М- количество исполь- 
зуемых характеристик), где каждая размерность может учитываться с некоторым 
весовым коэффициентом. Непосредственно сегментация осуществляется одним из 
распространенных алгоритмов сегментации (например, К-внутригрупповых средних 
или ОМА). Схема алгоритма представлена на рис. 6. 


ИА 
Матрица характеристик НН 


И 
Матрица характеристик? -—^__ Сегментация И 
№ мерных 
данных 
юнш 
ИМ 


Матрица характеристик М № 


Рисунок 6 — Схема алгоритма совместной сегментации 
(И’, — весовой коэффициент) 


528 «Искусственный интеллект» 372009 


Метод выделения областей на изображениях... лог 


Результаты экспериментов с алгоритмом сегментации 


Целью проведения экспериментов являлась проверка работоспособности 
разработанного алгоритма на имеющихся исходных данных, а также оценка по- 
грешности его работы. 

В качестве исходных данных для эксперимента использованы снимки участка 
картофельного поля, выполненные с разных высот со 2 по 6 июля 2007 года. На рис. 1 
приведены примеры исходных изображений, выполненных с высоты 100, 50, 15 
и 5 метров. 

Пример результата работы алгоритма приведен на рис. 7 (исходный снимок, 
выполненный с высоты 15 метров, приведен на рис. 1). 


Рисунок 7 — Пример результата применения разработанного алгоритма 


Полученный результат сегментации позволяет в автоматическом режиме выделить 
участки, на которых происходит развитие заболевания. Знание расположения таких 
участков позволит определить потребность тех или иных участков сельскохозяйственных 
полей в удобрениях и иных химикатах. А это позволит сделать сельскохозяйственные 
работы более эффективными и менее дорогостоящими. 

При решении построения карт особых областей возможно появление на карте 
участков, ошибочно отнесенных к особым областям. Это связано в частности с раз- 
мытием границ между частями исходного изображения (например, возникающий пе- 
реход цвета на границе изображений отдельных растений и почвы), которое может 
появляться, например, в результате ЛРЕС-сжатия исходных данных. 


Таблица 2 — Оценки погрешностей для данных различного пространственного 
разрешения 


Высота съемки, м 100 50 15 5 
№ эксперимента Численное значение оценки погрешности, % 
1 44,7 6,6 24,8 2,2 
2 4,6 26,6 4,6 3,9 
3 31,9 5,4 27 5,1 
Среднее 21,2 10,2 9,5 34 
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Указанная погрешность оценивается как относительная площадь ошибочных 
участков, полученная на исходном изображении, не содержащем особых областей, и 
зависит от пространственного разрешения исходных данных (оценки для различного 
разрешения приведены в табл. 2). 


Практическое применение алгоритма 


На основе разработанного алгоритма была предложена система принятия решения 
для программно-аппаратного комплекса внесения минеральных удобрений и других 
химикатов на сельскохозяйственные поля. 

Схема предложенного комплекса приведена на рис. 8. 

При работе по предложенной схеме внесение химикатов осуществляется следу- 
ющим образом: 

1. С использованием предложенного алгоритма вычисляются карты специальных 
областей (например, участки с развивающимся заболеванием растений). 

2. Построенные карты получают географическую привязку и сохраняются в 
базе ГИС для дальнейшего использования. 

3. Полученные карты используются при принятии решения о необходимости вне- 
сения того или иного количества удобрений на том либо ином участке сельхозугодий. 

4. Система управления внесением химикатов на основании имеющихся карт и 
данных, полученных в реальном времени, контролирует количество химикатов, вно- 
симых в почву, и направляет соответствующую команду системе внесения химикатов. 


Исходные Система принятия решений 
изображения 


Система 
Построение управления 
карт особых внесением 


областей химикатов 


Система 
внесения 
химикатов 


Предобработка 


Принятие 
решения о 


ЕЕ. внесении 

Пользователь химикатов 

Рисунок 8 — Схема программно-аппаратного комплекса для внесения минеральных 
удобрений и других химикатов на сельскохозяйственные поля 


В качестве данных, получаемых в реальном времени, могут использоваться: 
— данные глобальной навигационной спутниковой системы. В этом случае система 
управления, определяя с помощью навигационной системы, на каком участке поля 
она находится, вычисляет необходимое количество химикатов, исходя из данных кар- 
ты особых областей; 
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— данные от цветной камеры видимого диапазона. В этом случае система управления 
может в реальном времени корректировать данные карты особых областей, при этом 
принимая более точные решения, что увеличивает эффективность решения задачи. 


Заключение 


Полученные результаты по выделению фрактальных и текстурных признаков 
изображений сельскохозяйственных полей позволяют получить дополнительную инфор- 
мацию о состоянии растительного покрова. При этом результат в меньшей степени за- 
висит от условий освещения и наличия в кадре посторонних предметов, чем при исполь- 
зовании только цветовых признаков. Однако использование только этих признаков 
не позволяет достаточно четко выделить пораженные заболеванием участки. Для 
решения этой проблемы вводится алгоритм совместной сегментации, который для 
выделения особых областей использует также и цветовые характеристики изобра- 
жений для получения более точного результата. На основании этого метода и стро- 
ится программная система принятия решений о внесении минеральных удобрений и 
других химикатов на сельскохозяйственные поля. 
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В. Ганченко, О. Дудкт, А. Петровський, Б. Собковяк 

Метод видлення галузей на зображеннях з використанням фрактальних 1 текстурних характеристик 
знымюв вищого дозволу 

У стати пропонуеться метод обробки зображень, заснований на текстурних 1 фрактальних характеристиках 
зображень. Наведен! результати сум1сно! сегментацй з використанням розробленого методу при обробщ 
аерофотознмкв с1льськогосподарських полв та схема програмно-апаратно! системи прийняття рипення, що 
використовуе розроблений алгоритм обробки. 


Г. Сапсйенко, А. Роиайт, А. Регоу5Ку, В. зоБКоилаЕ 

Меод оЁ Агеа5 Эеесвоп оп Ппасе$ Ото Егаса! апд Техеига! СБагасе115с$ оЁ Ной Везош@йоп Ппаре5 
ТБе плаее ргосеззте тео фазе оп {ехга| ап4 Васа! сБагасет$Нс$ оЁ ппасе$ 15 ргорозе ш \е рарег. 
ТБе геза5 оЁ ]ошё зестещайоп изше фе 4еуе]оре тейо4 аге детопзиае ай ргосеззше аепа| 
рвоюзтарВ$ оРагиси ага] Не!4$. ТБе 4еуореа ргосеззте а1еогИбт 1$ аррНе4 ш 4ес151оп-таКте зузет, 
УСН 15 а раг оЁ зой\уаге-Баг\маге сотр[ех Юг азмсиЦаге аррИсайоп. 


Статья поступила в редакцию 09.06.2009. 
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